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Abstract 

The Hawking radiation in the vacuum of the spontaneous symmetry breaking in the 
gauge-Higgs- Yukawa theory is investigated by a general relativistic formulation of the 
the ballistic model. The restoration of the symmetry on the horizon and the formation 
of the spherical domain wall around the black hole are shown even if the Hawking 
temperature is lower than the critical temperature of the phase transition in the gauge- 
Higgs- Yukawa theory. When the Hawking temperature is much lower than the critical 
temperature, the domain wall closely near the horizon is formed. The wall is formed by 
the field dynamics rather than the thermal phase transition. 
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1 INTRODUCTION 



One of the most interestings in the quantum field theory with the general relativity is the 
thermal radiation from black holes, which is known as the Hawking radiation |2]. The 
temperature and the intensity of the radiation from the Schwarzschild black hole increase 
explosively at the final stage of the evaporation because the black hole loses its mass by 
the radiation and its temperature is inversely proportional to the mass 2J. Therefore the 
particle physics and the field dynamics including phase transitions around the black hole 
are especially interesting. Several authors discussed that the heating-up by the Hawking 
radiation can thermalize the neighborhood of the black hole and the thermal phase transition 
around the black hole can arise. Cline considered the QCD phase transition around the black 
hole whose Hawking temperature is greater than the critical temperature of the quark-gluon- 
plasma (QGP) transition He pointed out that some of the gamma ray bursts can be 
explained by the radiation from the QGP fireballs produced by the Hawking radiation from 
the primordial black holes. We showed the thermal electroweak (EW) phase transition around 
the black hole and the formation of the spherical domain wall which separates the symmetric- 
phase-region from the broken-phase background E] • We discussed several mechanisms of 
the baryon-number-production by the spherical domain wall around a black hole and proposed 
the cosmological baryogenesis scenario by the primordial black holes [U El E] • 

The thermal phase transition on some local volume requires the local thermal equilib- 
rium which is confirmed by the enough interaction rate in the considered volume. The mean 
free paths of the radiated particles are much longer than the Schwarzschild radius. If the 
thermal phase transition is caused by heating-up by the Hawking radiation, the radius of 
the wall should be much larger than the Schwarzschild radius. Hence the formation of the 
wall by the thermal phase transition needs much higher Hawking-temperature than the crit- 
ical temperature of the transition [HEl- Here, it is natural to ask if the Hawking radiation 
whose temperature is similar to the critical temperature of the transition infiuences the Higgs 
vacuum-expectation- value (vev) around the black hole. In our previous work [HIIH], we con- 
sidered the Hawking radiation with such a temperature in the gauge-Higgs- Yukawa theory, 
e.g., the EW theory or the Grand Unified Theory (GUT). We proposed the ballistic model as 
an effective description of the system. By using the ballistic model, the formation of the wall- 
structure of the Higgs vev was shown when the Hawking temperature is equal to or greater 
than the energy scale of the gauge-Higgs- Yukawa theory. The wall-structure is formed by the 
field dynamics rather than the thermal phase transition. 

The ballistic model is described by both the action for the Higgs scalar field and that for the 
relativistic point particles as the Hawking-radiated particles from the black hole. The ballistic 
description of the radiated particle is valid because the mean free paths determined by gauge 
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Figure 1: The balance of pressures — our basic idea for the wall-formation in the ballistic 
model. The parameter r is the distance from the center of the black hole and | {4>{r)) \ means 
the absolute value of the Higgs vev. Once a wall-structure of Higgs vev (thick curve) is 
formed, mass of the particle is depending on the distance r. A part of the momentum of the 
Hawking-radiated particle Ap = p{r) — p{r + Ar) is transfered to the wall-structure because 
the momentum is depending on the distance r in the wall-structure: p{r) = ^ E"^ — m{ry. 
The transfered momentum causes a radiation-pressure acting on the wall. A local balance 
among the radiation-pressure, the wall-tension and the potential- force can stabilize the wall- 
structure. The dotted thick curve means the wall-structure near the horizon. The radiation 
pressure near the horizon is enhanced due to a general relativistic effect and forms the wall- 
structure near the horizon. 
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interactions are much longer than the Schwarzschild radius [H |H] and the mean wavelength 
of the radiated particles approximately given by the Schwarzschild radius is shorter than the 
length-scale of the wall-structure [3 EI- The basic idea of the wall formation in the ballistic 
model is schematically shown in Figure Q When once a radius-depending Higgs-vev-structure 
I {4>{r)) I, namely the wall-structure, has been formed, the Hawking-radiated particles run rush 
up the slope on the wall-structure and they push the wall-structure outside. Therefore there 
arises a pressure, named the Hawking-radiation- pressure, which is acting on the spherical 
wall-structure to expand its radius. On the other hand the wall-structure also feels its wall- 
tension and the Higgs-potential-force, which are acting on the wall-structure to shrink its 
radius. The balance among these pressures keeps the wall-structure stationary. 

In the previous work jH] , the black hole is assumed as a simple source of the particles with 
thermal spectrum and general relativistic (GR) effects are omitted. In this paper we fully 
consider the GR effects in the ballistic model and show the formation of the wall-structure 
if the Hawking temperature is not only greater than the energy-scale of the gauge-Higgs- 
Yukawa theory but also smaller than the energy-scale. We assume that (i) the motion of 
the particles radiated from the horizon obeys the geodesic on the Schwarzschild space-time 
and (ii) the observer at the infinite distance detects the disk-image with uniform intensity 
as an image of the radiated particles. By constructing the general relativistic formulation 
of the ballistic model we derive the effective Higgs potential including full GR effects, which 
determines the Higgs vev structure around the black hole. The effective potential is depending 
on the differential particle-density-distribution dE x //{E, r) of both the particle-position r 
and the particle-energy E. The analysis of the particle-trajectories radiated from the horizon 
gives us the density-distribution. The density-distribution dE x Af{E, r) is depending on the 
differential particle-flux dE x duj sin uj x f[E,uj) on the horizon because the initial conditions 
for the particles are given by the flux f{E,uj) on the horizon. The parameter is a zenith- 
angle of the particle-radiation on the horizon. The elevation-angle of the radiation becomes 
(7r/2 — uj). The particle-fiux for high elevation-angles is derived by the two assumptions 
above-mentioned because particles radiated with the high elevation-angle can be detected by 
the observer and their trajectories can be traced. On the other hand, the particle-fiux for 
low elevation-angles cannot be determined by the assumptions because the particles radiated 
with a low elevation-angle return into the horizon and are not detected by the observer. 
By the analytic continuation of the particle-fiux from the high elevation-angle to the low 
elevation- angle, we can obtain the particle flux for all angle. Then we can derive the effective 
potential and the effective differential equation for the Higgs scalar vev around the black hole. 
The suitable boundary conditions for the differential equation are as follows; (i) the Higgs 
vev approaches the ordinary Higgs vev (0) = f/-\/2 which minimizes the bare Higgs 

potential V{(f)) when the distance from the black hole approaches infinity r ^ oo and (ii) 
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the Higgs vev should be finite on the horizon. The field equation is numerically solved with 
the boundary conditions and the wall-structure of the Higgs vev around the black hole is 
obtained. 

Our resultant wall-structure has two interesting properties as follows. The Higgs wall 
structure is formed even if the Hawking temperature is much smaller than the energy scale of 
the gauge-Higgs- Yukawa theory. When we consider such a low temperature black hole, there 
arises a wall-structure closely near the horizon. The Higgs vev on the horizon becomes zero, 
namely, the symmetry broken-down spontaneously by the Higgs potential is restored on the 
horizon. These properties are not found in the previous work [S^ and are essentially caused 
by the GR effects on the Hawking-radiated particles. 

In Section |21 the general relativistic formulation of the ballistic model is given and the 
general form of the effective Higgs potential Ves is obtained. The effective potential Ves 
is depending on the differential density-distribution dE x Af{E, r) of the radiated particle. 
In Section El various patterns of the particle-trajectories are analyzed and the form of the 
density-distribution dE x Af{E, r) is obtained as a functional of the differential flux on the 
horizon dE x duosmuj x f[E,uj). In Section H] the flux dE x duj sinuj x f[E,uj) is discussed. 
In Section El the concrete form of the density-distribution dE x f/[E,r) is obtained by using 
the results of Section El and Section HI In Section El the differential equation for the Higgs vev 
and the boundary conditions are discussed and the formation of the wall-structure is shown. 
In section [71 we provide a conclusion and discussions. In Appendix A another argument of the 
flux on the horizon dE x dujsinuj x f[E,uj) is given. In Appendix B the energy-density of 
the Hawking-radiated particles around the black hole is calculated. In Appendix C analytic 
solutions near the horizon are discussed. 

2 Ballistic Model with General Relativity 

The ballistic model is an effective theory for computing the structure of the Higgs vacuum 
expectation value (vev) around a black hole radiating particles in the vacuum of the gauge- 
Higgs- Yukawa theory [Sj. The Hawking-radiated particle can be regard as ballistic because 
the radiated particle has a longer mean free path and has a shorter wavelength than the 
radius of the wall-structure which we will discuss later. The ballistic particle is identifled 
with a relativistic point particle with interactions. In the ballistic model the Higgs scalar 
fleld is separated into the Higgs vev and the ballistic Higgs particles [Sj. The Higgs vev is 
depending on the position but is not depending on the time. The ballistic Higgs particle 
describes the Higgs propagating mode in the background of the Higgs vev. The ballistic 
model consists of the time-independent Higgs scalar vev (f){x) and a set of the trajectories of 
the ballistic particles (s)}, where i denotes an index of each particle and s is a parameter of 
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the trajectory. The set of the balhstic particles {l/f (s)} contains all radiated particles around 
the black hole, which also includes the ballistic Higgs particles. The fundamental action for 
the ballistic model is given by a combination of an action for the Higgs scalar vev (j){x) with a 
Higgs potential V{(f)) and an action for the trajectories {yf (s)} of relativistic point particles 
as the ballistic particles: 



S[<l>,y] 



J2 / dsYi\(f){yi)\Jg^^yl'y';; 



1) 



where we have defined := dy^/ds and the summation is performed over all ballistic parti 
cles. The background space-time is given by the Schwarzschild metric: 

ds'^ = +F{r)dt^ - F-\r)dr'^ - r^dO^ - r^sin^^rfy?^ 

with the Schwarzschild factor: 



F(r) : = 

r 

The bare Higgs potential in the action is given by the double-well form: 



(2) 
(3) 



.2 J 2 



+ 



(4) 



which has a minimum at |<^| = v/V2 and the constant /x^ > is the Higgs mass, therefore, 
the Higgs vev without the Hawking radiation is given by | (0) | = v/\^. We have defined an 
extended- Yukawa-coupling-constant Yi for the point particle i. The mass of the particle i is 
proportional to the Higgs vev and is given by 



rrii 



^.1(0) 



(5) 



where the mass relation in the gauge-Higgs- Yukawa theory is reproduced. Especially mass of 
the particle i in the vacuum is given by rrii = Yi\ {(p) \ = Yiv/^/2. The ballistic Higgs particles 
does not obey the relation (0), however, by choosing 

iHiggs = VS^ (6) 

as the extended- Yukawa-coupling-constant for the ballistic Higgs particles, the action ^ 
results in the correct effective potential (see Appendix in [H]). The formulation presented 
above is a simple general relativistic extension of the original ballistic model in the previous 
work [H]. 

By a gauge fixing ds = dt = y\ and employing the Schwarzschild coordinate (t, r, 6, (p), the 
action becomes 



S[4>,y] 



X 



dtdrdOdip sin 6 

Y,m,r,e,^)\ 



r-\de(j)y - r- 



sm 



1 5{r-y-)5{d-y')5{^-y^) 



sin 6 



(7) 
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where we have defined a general relativistic gamma factor for particle yi as 



F - F-'iij^' - r'iy^ - sin e{ytf] . (8) 



We assume spherical symmetry of the Higgs vev 0(r) due to the spherical symmetry of 
the metric and the Hawking radiation. The equation of motion for the Higgs vev with the 
condition dtcj) = becomes 

where we have defined a Laplacian with spherical symmetry 

A0 := + ^dr{r^Fdr(p) . (10) 

Because the metric and the Higgs vev as a background of the ballistic particles have an 
invariance of the time-evolution and have spherical symmetry, there are two constants of 
motion for each ballistic particle, namely, the energy Ei and the angular momentum L,. 
By using spherical symmetry we select the coordinate system individually for each particle 
trajectory to keep = 0. The motion of each particle is described on each plane (?/[, yf ). By 
using constants of the motion, the equation of motion for the ballistic particles becomes 

= Ymyl)\l^{t)F{yl), (11) 
= r,|0(y[)|7.(t)r2y^. (12) 

By substituting (fTTj) into © and using property of the Dirac delta function, we obtain 
. , _ ^dV 19101^,,, ^,J<Pir)\ S{r-yr)6{e-y^)6{^-yf) 

^'^ " 2 90 + ' ^^"^^"B^ ^^^e ■ ^^^^ 

We define the effective potential as 

VM*,y) := VW + ^/f(r)E>^'^ ''''"^-'''!:f (14) 

then the equation of motion for the Higgs vev becomes 

A0 = (15) 
^ 2 90 ^ ' 

Instead of considering the effect of each ballistic particle, we will statistically treat the 
group of the particle of the same kind. We adopt a differential distribution of particle- 
number-density dE X Mf{E,r) for particle species / with energy E at the position r. The 
density-distribution should has the spherical symmetry because of the absence of a special 
direction of the Hawking radiation. The effective potential (jl4j) becomes 

1 r d W 

Fefr(0,r) = y(0) + -02F(r)^r/y — Ar^(i?,r), (16) 
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where the summation in the effective potential is performed over all particle-species {/} 
in the gauge-Higgs- Yukawa theory and Yf is the extended- Yukawa-coupling-constant for a 
particle-species /. We rewrite the effective potential (fT^ as 

Kff(0,r;Ar) = +l;,2^(r;A^)02 + (17) 
where we have defined the effective Higgs mass as 

t^Ur-J^) = V + nr)Y.yf j -Y^fiE,r). (18) 

The effective Higgs mass ()18p which governs the effective potential ()17|) is a functional of the 
density-distribution A//(-E',r). Therefore we have found the Higgs vev structure around the 
black hole is determined by the density-distribution of the Hawking-radiated particles. 



3 Density of the Ballistic Particles I 

In this section we calculate the density-distribution A//(-E', r) of the ballistic particles radiated 
from the black hole. To determine the density-distribution we assume as follows: 

(i) The particles are radiated from the horizon and obey the equations of the motion in 
(fTT|) and (fT^ on the Schwarzschild space-time Q. 

The configuration of our setup is schematically shown in Figure El A part of radiated particles 
reach infinite distance. The initial condition for the particle is given by the position {6, (p) 
where the particle radiate on the horizon, the radiation-angle {uj, ip) on the position and the 
energy E which is defined in ()11|) . is a zenith-angle defined for < < 7r/2 and ip is a. 
azimuthal angle defined for {] < ip < 2t[. The zenith- angle uj is defined by the angle which 
the trajectory of the particle and the perpendicular of the horizon make on the horizon. The 
zenith-angle uo is defined on the Schwarzschild coordinate. The elevation- angle of the radiation 
becomes (vr/2 — uj). 

We put the differential number-flux of the Hawking-radiated particles per an unit area on 
the horizon with the Schwarzschild coordinate 

dTf = f f{E,uj) X dE X duj sin UJ dip (19) 

which depends on the energy E of the particle and on the zenith-angle u of the radiation. 
The flux (fT^ is independent of the position {6, 0) and of the azimuthal angle ip due to the 
spherical symmetry. The form of the particle-trajectory is determined by the particle-energy 
E and zenith-angle u of the radiation as the initial conditions, and is also determined by the 
mass of the particle, i.e., the product of the extended- Yukawa-constant Yf and the Higgs vev 
|0(r)| as the back ground. Then the density-distribution Aff{E, r) is obtained by the analysis 
of the particle-trajectories and is depending on the flux ff{E,uj). 
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Figure 2: The configuration of our setup — tlie coordinate {9, 0) of the horizon, the radiation- 
angle {oo,il)) on the horizon and the screen at the infinite distance, uj is the zenith-angle and 
%l) is the azimuthal-angle of the radiation on the horizon. The elevation-angle of the radiation 
becomes ('iT/2 — uj). The screen is parameterized by the polar coordinate (Z^,x), where D is 
the distance from the center of the screen and x is the angle. D is identified as the impact- 
parameter of the particle-geodesic. The observer at the infinite distance finds an image of a 
disk as a detection of the Hawking radiation. Dc is the radius of the disk-image, which is 
depending on the ratio of the particle-energy E to the particle-mass m. In this figure the 
0-axis and x-axis correspond, however, this correspondence is dispensable. The strength of 
the Hawking radiation is independent of {9, (p) and of ip due to the spherical symmetry of 
the black hole and of the Hawking radiation. The disk-image is independent of x because 
of the spherical symmetry. Each value of -D (< Dc) is one-to-one corresponding to lj by the 
particle-geodesic when the particle-energy E (> m) is given. 
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3.1 Classification of Particle Trajectory 

We rewrite the motion-equations and ()12|) for the balhstic particle into the differential 
equations: 



EJ {y^-f \ E 

We will abbreviate the particle-index i or flavor-index / in equations for simplicity. By 
combining the equations (j2(Jj) and (j2H) a relation for the zenith-angle uj is obtained as 



dy^ 

tana; = Tbh 



dy- 



E tbh 



LI 

(22) 



By using the relation ()22j] we employ the zenith-angle a; as a constant of the particle-motion 
instead of the angular momentum L. We approximate by 

F|0(r)| ~ m (23) 

in the motion-equation ()2()j] for simplifying the analysis, where m is the mass of the particle 
in the ordinary vacuum = v/ \/?) . By this approximation the backreaction from the Higgs- 
vev-structure into the particle-trajectories is ignored and the only influence of the radiated 
particles over the Higgs vev is considered. The equations of motion for the ballistic particle 
in (I2II and (jHI) becomes 



{y'f = -WM, (24) 

rBH 



yt> = tano^^F, (25) 



where we have defined an effective potential for the particle-motion 



tbh \^ ^ fm 



WMV) ■■= -F' {1 - tan^c {^j F - [-^fF}. (26) 

The effective potential is parameterized by two constants {m/E,uj). 

According to the function-form of the effective potential ()26|1 . the parameter space {m/E, u) 

is divided into the six regions; the Region-I, the Region-IIa, the Region-IIb, the Region-IIb+, 

the Region-IIb— and the Region-IIc. The division of the parameter space is shown in Figure 

ini On each region the effective potential Wes has different forms displayed in Figure |31 In 

the Region-I the particle radiated from the horizon runs away from the horizon to the infinite 

^This approximation is available because the form of the final results (j>{r) in Figure |51 are approximately 
flat. 
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distance because the effective potential Wes never cross tfie zero level except for the horizon 
f = ^BK- In the Region-II's the radiated particle always turns its r-direction at each point 
satisfying Wesir) = and returns into the horizon r = tbh- The trajectory of the radiated 
particle is connected with the horizon and is described by each thick curve in Figure IH In the 
Region-IIa and in the Region-IIc there is a trajectory disconnected from the horizon (the dot- 
ted curve in Figure |3)). The trajectory disconnected from the horizon in the Region-IIc means 
the revolution-orbit around the black hole and that in the Region-IIa means the gravitational- 
scattering-trajectory by the black hole. Such trajectories disconnected from the horizon are 
not relevant to our subject. 

In the Region-II's the relation amang the turning point rturn, the radiation angle u and 
the ratio m/E is given by Wesirtum) = and it can be rewritten as 



tan^ uj 



tan Utv 



m 

'e 



(27) 



where we have defined a function 



tan Utu 



m 



^BH 



(28) 



We display the form of the function tantutum in Figure El When < m/E < 1, the relation 
is valid for tbh < '"tu 



< rciim/E) as the Region-IIa. When 1 < m/E < y9/8, the 
relation (^Tj) is valid for tbh < '"turn < rdim/E) as the Region-IIc and the Region-IIb-|- and 
for Tcs < Tturn < rc2{rn/E) as the Region-IIb— . The relation (P7j) is not valid for Vd < rturn < 
fcs{fn/E) because the relation indicates the turning radius of the revolution-orbit. We note 



that any particle with 1 < m/E < J9/8 does not turn at the radius Vd < r < Vcsim/E). 



When y9/8 < m/E, the relation is valid for tbh < ^tum < rc2{m/E) as the Region-IIb. We 
have defined the first maximum turning radius 



m 



'^BH 



/ m V 



m\2 3 / 8 /m 
e) +4V^-9b 



(29) 



for < m/ E < J 9/8 and its dual radius 



/ m\ 

'"clDual ( ^ ) 



'^BH 



l-(f)' 



4 \e) 4V^ 9\e) 



my 



(30) 



for 1 < m/E < J 9/8. These are given by the solutions of the quadratic equation 



d 



Or, 



■ tan LJt^ 



0. 



(31) 



turn 



'^turn — "^c 
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We have also defined the second maximum turning radius 



rc2 



m 



( HL 



^' (f) -1 

for 1 < mjE as a solution of tan^a;(rc2, wi/i?) = and have defined the splitting radius as 



(32) 



E 



'my 
9 \E) 



(33) 



for 1 < mjE < y9/8 as another solution of tan^ ciJturn(^cs) = tan^ co'turn('^ci)- The critical 
curve u = udm/E) which separates the Region-I from the Region-IIa in FigureElis given by 



tan uJr 



2 I ^\ 

tan tUturn ( Tcl, — 1 



7 imV^ l-l + 27[l-|(f) 



3/2 



l-(f 



(34) 



where the right hand side is the function defined in (j28j). This critical curve also indicates 
the boundary which separates the Region-IIc from the Region-IIb— . This critical curve is 



defined for < m/E < ^9/8. The other boundary which separates the Region-IIc from the 
Region-IIb+ is given by 



tan UJcBnal 



2 / rn\ 

tan UJ [ TclDual, -gT ) 



7 
2 



m\ 2 1 
e) +8- 



27 



1 3/2 



l-(f 



(35) 



which is defined for 1 < m/E < j9/8. 



3.2 Density Distribution as Functional of Flux 

We consider the density-distribution J\ff{E, r). The density-distribution is a functional of the 
flux f{E, uj) on the horizon and is also a function of both the parameter mf/E and the position 
r, then we rewrite the density-distribution as //{m/ E,r] /) in this subsection. The particle- 
density at the position r is the superposition of the density of the particles which can reach 
the position r. When we put values of both the radiation angle u and the parameter m/E, 
we find whether the particle can reach the position r or not. A contribution to the density by 
a particle is inversely proportional to its velocity y^. Therefore the particle- number-density is 
calculated by integrating out the product of the flux f{E,uj) on the horizon, the geometrical 
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Figure 3: Parameter space {m/E,uj) of the effective potential PVefr for the balhstic particle. 
The space is classified into six regions; Region-I, Region-IIa, Region-IIb, Region- 1 lb -|-, Region- 
lib— and Region-IIc. The typical forms of the effective potential corresponding to the regions 
are shown in Figure El The radiated particle runs away from the horizon in the Region-I. In 
the other regions (Region-II's), any radiated particle returns into the horizon. The Region-I 
is bounded on both the curve uJc{m/E) in (jMj) by the Region-IIa and the line m/E = 1 by 



the Region-IIb—. We note tanc<Jc(0) = \J21/A and tanco'c(l) = 2. The Region-IIc is bounded 
on the curve cJcDuail^^Z-^) in ffH^ by the Region-IIb+, on the curve uJcim/E) by the Region- 
IIb— and on the line m/E = 1 by the Region-IIa. The Region-IIb is bounded on the line 

m/E = 
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Region-I Region-IIb- 



Figure 4: Typical forms of the effective potential Wes{r) in ()26p for the six parameter regions. 
The figures are arranged for corresponding to the position of each region in the parameter 
space in Figure El The curve is the function- form of the effective potential Wes{r) for the 
ballistic particle and the motion of the particle is allowed for Wes{r) < 0. The thick curve 
means the trajectory of particle-motion which is connected to the horizon. The motion 
of Hawking-radiated particle is described by the thick curves. The dotted curves mean the 
particle-trajectories disconnected to the horizon, i.e., the gravitational bending in the Region- 
Ila and the revolution-orbit in the Region-IIc, and is not relevant to our subject. In the 
Region-I the radiated particle runs away from the horizon into the infinite distance. In the 
Region-II's the radiated particle turns its r-direction at the radius r = rtum and returns into 
the horizon r = Tbh- 
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Figure 5: The function-forms of tan tuturn ('"turn, '^/-E) defined in for various values of 
m/E > 0. The forms of the curve are classified into three classes < m/E < 1, 1 < m/E < 
^9/8 and y^9/8 < m/E. The curve gives us the relation between the turning radius rtum and 
the radiation- zenith-angle tanu; = ta.nujt^j.^{rt^j.^,m/ E). For given u, the solution r nearest 
the horizon becomes the true turning radius rtum- Then solid curve is the turning radius 
rturn and the dotted curve means a fake. For < m/E < ^9/8, the solid curve is bounded 
on rturn = fci which is the first maximum turning radius whose value is given in ()29|). On 
f^turn = 0, the solid curve is also bounded on rtum = which is the second maximum turning 
radius (j32j) . For 1 < m/E < \j9/8, we have another bound on rtum = '"cs which is the splitting 
radius 
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factor [r^Yi/i')'^ ^^e inverse velocity l/y^' with the integral variable u. Here we define the 
function 

KtimlE,r,l) :^ (^^)' jy^... ^^A^ ,36) 

for convenience. Because a part of the parameter space {m/E,uj) contributes to the density- 
distribution J\f{m/ E f) at the position r, the integral interval [^71,072] for the density- 
distribution Af{m/E, r; /) should be carefully chosen, i.e., the integral interval [ui, UJ2] depends 
on the position r and on the parameter m/E. The density-distribution is given by the sum 
of contributions from the parameter regions: 

A/'(|,r;/) = M(|,r;/)+Mia(|,r;/)+Mib(|,r;/)+Mib±,c(|,r;/)(37) 

Afi{m/E,r; f) is the density-distribution contributed by the Region-I, Afua{m/ E,r; f) is that 
by the Region-IIa, Afiih{m/ E,r; f) is that by the Region-IIb and Niih±,c{jn/ E^r] f) is that 
by the combination of the Region-IIb—, the Region-IIb+ and the Region-IIc. These density- 
distributions are given by the function with a suitable integral interval [co'i,ti;2]. The 
density-distributions Afi{m/E,r; /), Afiu{m/E,r; /), J\fuh{m/E,r] f) and A/iib±,c("^/-E, r; /) 
are considered as follows respectively. 

1) The density-distribution contributed by the Region-I. 

This parameter-region is given by < m/E < 1 and < u < uJc{m/E) (see FigureEl). The 
radiated particle in this region runs away from the horizon r = tbh to the infinite distance 
r = 00 and never turns its r-direction (see Figure llj). Therefore the integral interval is given 
by < < uJc{m/E) and this interval does not depend on the position r. The density- 
distribution becomes 

^-^(f ) ^. f\ (0 < f < 1 and TBH < r < 00) 

^ . (38) 

(others) 

2) The density-distribution contributed by the Region-IIa. 

This region is given byO < m/E < 1 and uJc{m/E) < uo < tx/2 (see Figure Ej). The 
radiated particle runs away from the horizon r = tbh; turns its r-direction at radius r = rtum 
which depends on the zenith-angle uj and returns into the horizon, (see Figure E}. Then 
the integral interval [co'i,co'2] depends on the position r. The particle within the parameter 
region uJc{m,/E) < uj < uj^^„^{r, m/E) can reach the position r and contributes to the density- 
distribution Afiiaini/E, r; /) because a particle with the zenith-angle u = ti;turn('", m/E) turns 
its r-direction at the position r. On the other hand, the particle within the parameter region 
|^turn('"5 m/E) < UJ <Tx /2 cauuot reach the position r and does not contribute to the density- 
distribution M\\a.{m/ E f). Therefore the integral interval should be uJcijn/E) < uj < 
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^tnrnij','^/ E) which is depending on the position r and the density-distribution becomes 
Mia(|,r;/) = 



<T-) " i' ^' ^ (O < f < 1 and TBH < r < r,, (f )) 



(others) 

We need the factor two in this form because both the go and the return of the particle 
contribute to the density-distribution. Any particle in the Region-IIa does not contribute to 
the density-distribution for r > rd because the maximum turning radius in the Region- 
IIa is given by rd and particle does not turn on the radius r > r^i (^)- Then the 
integral form in (IHU)) is restricted to Tbh <r Kr^i (^)- For example we obtain tantUc = \J^- 
and = (3/2)rBH when we consider massless particle m = 0. We obtain tancUc = 2 and 

= 2rBH for the particle with m/E = 1. 

3) The density-distribution contributed by the Region-IIb. 



This region is given by ^9/8 < m/E and < c<j < 7r/2. The particle whose zenith-angle 
satisfies < u < cjtum (r, ^) can reach the position r because the radiated particle with the 
zenith-angle u = ootum (r, ^) turns its r-direction at the position r. The density at r is given 
by the integration over the radiation angle < uj < ciJturn (j', We obtain 

, I 2 X K--^^-^ (|,r; /) (4 < f and tbh < r < r,, (f )) 

A/'iib(^,r;/) = ^ (40) 



E- 



(others) 



The integral form in pn)) is restricted to tbh < r < (j^^ because the maximum turning 
radius in this region is given by rather than rd (^^^ (see the curves with ^9/8 < m/E 

in Figure El)- For example, we have = 9rBH for the particle with m/E = ^9/8 and we 
have — > 1 for the particle with the limit m/E —>■ oo. 

4) The density-distribution contributed by the Region-IIb± and the Region-IIc. 



The parameter-region is simply given by 1 < m/E < ^9/8 and < uj < 7r/2, however, 
the analysis becomes a little complicated and more careful treatment is required because the 
region is composed of three Regions (IIb+, lib— and lie) due to the revolution-orbit in the 
Region-IIc, which is irrelevant to our subject. Each particle radiated from the horizon with 
zenith-angle u turns its r-direction at the position rtum and returns into the horizon. The 
relation between u and rtum is given by the equation ()27|) (also see Figure EI). In the region- 
lib- (u < u;c(f ) := t^turn('^ci, f )) and in the region-IIb+ > a;cDuai(f ) := ^turn('^ciDuai, f )), 
we have a unique solution rtum of the equation (|27p for the fixed u. On the other hand, in the 
region-IIc < uj < tx;tum('"ciDuai)) there are three solutions of the equation (j77|l (see Figure 
0] and Figure 13)) because the revolution-orbit around the black hole is also allowed. The 
solution nearest to the horizon should be selected as rtum because the revolution-orbit which 
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is disconnected with the horizon is not relevant to the Hawking radiation and the trajectory 
of the Hawking radiation is connected to the horizon. 

When the zenith-angle uj becomes larger from zero to uJi.{m/ E) in the considering region 



1 < m/E < ^9/8, the turning radius rtum becomes smaller. When uj equals to oOcim/E), 
the turning radius rtum discontinuously jumps from the radius r^s to the smaller radius Vd 



(see the curve with 1 < m/E < ^9/8 in Figure Ej) because the parameter region is changed 
from the Region-IIb— into Region-IIc. At the moment the trajectory is separated into the 
two disconnected parts and one of the trajectories becomes the revolution-orbit which is 
irrelevant (see the difference between Wesir) in the Region- Hb— and that in Region-Hc in 
Figure|3)). When u becomes large from ujc{m/ E) to 7r/2, the turning radius rtum becomes small 
continuously and approaches the horizon radius tbh- Noting special happens in the transition 
from the Region-IIc into the Region-IIb+ except for the disappearance of the revolution-orbit 
which is not relevant. 

To consider the density-distribution at the position r, we should consider which particle 
can reach the position r. When rcs(^) < r < rc2(^), the particle belongs to the Region- 
IIb—. The particle radiated with the zenith-angle uj = cuturn turns its r-direction at 
the position r, then the particle with < uj < ujtum (j'^ can reach the position r. The 
density at r is given by the integration over the radiation angle < uj < cutum (r, ^) • When 
'^ci(^) < T < rcs(^), the particle still belongs to the Region-IIb—. The radiated particle 
does not turn in the region rci(^) < r < rcs{^) and always turns at r = rcs(^). Then 
the integral-interval should be < < ujtum (j'cs, The integral- interval can be rewritten 
as < < uJc{^) by using the definition of uj^i^) in (jHl)). When tbh < r < rci(^), the 
particle which belongs to the Region-IIc or the Region-IIb+ begins to contribute the density- 
distribution. The situation is essentially same to the Region-IIb— for r^si^) < r < rc2(^) 
and the integral interval is also given by < < uJtum ■ The density-distribution for 



A/iib±,c i-^^^'^f 



■ (41) 



1 < m/E < y9/8 is summarized as 

(rc2 (f ) < r < oo) 

2 X (f ,r; /) (r,3 (f ) < r < r^, (f )) 

2 X (f , r; /) (r,i (f ) < r < r.^ (f )) 

2xAC*-"('-'^^(f,r;/) (rBH<r<r,i(f)) 

The resultant form ()4H) has a little complicated conditions for r, however, the form is a 
continuous function of r because the integral interval continuously changes as r changes. 



18 



4 Radiation Angle Distribution 



In this section we consider the particle-flux on the horizon f{E,uj) in p9|) . which is required 
to evaluate the particle density-distribution Af{m/ E,r; f) in (jHTj). 

As a consequence of the assumption- (i) which is mentioned in the previous section and of 
the spherical symmetry, the observer at the infinite distance finds a disk image as a detection 
of the Hawking-radiated particles. In Figure |2l the correspondence of a particle-geodesic to 
a point in the image of the radiation on the screen at the infinite distance is schematically 
shown. The observed image is equal to the image on the screen. The screen is parameterized 
by the polar coordinate {D,x), where the parameter D is distance from the center of the 
screen. The parameter D can be regard the impact parameter of an absorption-process into 
the black hole because the absorption of the particle is regarded as the inverse process of the 
Hawking radiation. A relation 

§ = v^D, (42) 

is derived from conservation of the angular momentum of the particle, where velocity of the 
particle at the infinite distance 



is defined. The particle which reaches the infinite distance satisfies a condition E > m. By 
combining (|43p with (|22p we obtain the relation between the radiation- zenith-angle u and the 
impact parameter D as 

D(-,u;) = TBK . 44 

\E J Voo 

In order for a radiated particle to reach the infinite distance, the radiation- zenith-angle should 
satisfy < tancu < ta.nujc{^). This restriction for the radiation-zenith-angle u is correspond- 
ing to the restriction for the impact parameter < D < Dc{^). By using the relation between 
u and D in (j44|) . we obtain the radius of the apparent disk- image on the screen as 

dJ-) := TBH -. 45 

\EJ Voc 

Therefore the impact parameter D, namely the position of the image, is one-to-one corre- 
sponding to the radiation-zenith-angle u ii D < Dc{m/E) and E > m are satisfied. 
The apparent disk 



means the absorption-cross-section of the black hole for the particle in the ballistic picture 

27„ 



(eikonal limit). The area of the disk is depending on m/ E. For example we have Dc = y -^tbh 
and (Tbh = X'^'^bh ^'^^ massless particles. 
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Because the Hawking radiation for the observer at the infinite distance is regarded as an 
inverse process of the particle-absorption into the horizon, the differential luminosity of the 
Hawking radiation is given by the product of the absorption-cross-section and the thermal 
fiux from the black body per unit area: 

,^Bn('^)x^-fx-^fT,^iE)47iE'v^dE{E>m) 
dL{TB}i,m,E) = { \EJ 4 [2tiY , (47) 



{E <m) 



where 



fT^^{E) := (48) 

e^BH ± 1 

is the Fermi-Dirac or the Bose-Einstein distribution-function for the Hawking temperature Tbh 
and g is the degree of the freedom for the particle on the temperature Tbh- The absorption- 
cross-section (Tbh in the luminosity dL is depending on the particle-energy E. This 
dependency is referred as "the gray body factor" of the Hawking radiation mi2lE|- We have 
a relation between the fiux on the horizon f{E, uj) in (fT^ and the differential luminosity ()47|) : 

dT = 47rrBH / dusinuj / dxf{E,uo)dE = dL{TBH,m,E). (49) 
Jo Jo 

This relation gives us a restriction for the flux f{E,uj). 
To determine the flux f{E,uj), we assume 

(ii) the brightness of the disk-image for each particle-energy appears to be uniform over 
the disk-image for the observer at the inflnite distance and the energy-spectrum of the 
disk-image is given by the perfect black body radiation with the Hawking temperature. 

By the above assumption and the relation ()44|1 we have the relation between the differential 
brightness of the disk (jTTj) and the flux on the horizon as 



T (lo^^^««^^^ ^TrE^^dE X (dDDdx) 



4nrl^ X 



1 1 



4n cos^o; 4 {271^'^ 



]j^MAE)A^E'dE 



X (dtusintudx). (50) 



By using the spherical symmetry of the Hawking radiation, [du sin udx) in the right hand side 
of ()50j) can be rewritten to {duism-udip)"^. By comparing with the deflnition of the particle-flux 
in (|T^. we obtain the particle-flux on the horizon per unit area 



9_ 

IGvrcos^a; (27r) 



f{E,u;) = ——^Xj£-fT,^{E)4nE'. (51) 



^dx is equal to d4> due to the configuration in Figure|5] dcj) can be regard as dip because both the space-time 
and the Hawking-radiation are spherical symmetric then 0-axis is exchangeable for ^p-axis. 
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The resultant flux satisfies the restriction in 

We should note that the resultant flux ()5H) is defined for the particles which belong to the 
Region-I because the particles in Region-II's cannot reach the observer. Namely, the flux (j^Tjl 
is only defined for high elevation-angle^ < uj < Uc and high particle-energy E > rrif. We 
need the flux f{E,uj) defined for all elevation- angle < u < n /2 and for all particle-energy 
< < oo to calculate the particle number density N'{m/E,r] f). To solve the problem 
we perform a analytic continuation of the flux, namely, we extend the domain of the flux 
1)5111 to all elevation- angle < u < 7[/2 and all particle-energy < < oo with keeping 
the function form of the flux (jKT|) . The analytic continuation can be performed C°°-smoothly 
because no singularity exists in the new domain of definition. The resultant flux ()3T|) defined 
for < u; < 7r/2 diverges on uj = tt/2. 

The radiation-flux on the normal black body with temperature Tbh is given by 

fiE,u;) = '-^v'xj£yJr,^iE)4nE^ (52) 

Our resultant flux on the horizon (jHT|l is different from the flux (jK^ . Especially the dependence 
on the radiation- zenith- angle u is quite different. We find the low-elevation- angle-dominance 
of the flux on the horizon. The dependence on the particle- velocity is also different. 



5 Density of the Ballistic Particles II 

We concretely calculate the particle-density-distribution ()H7)) by using the resultant flux on 
the horizon We prepare several functions to calculate the density-distribution (jSTj). By 
substituting the flux f{E,u) resulted in into the functional J\f^^{m/ E,r; f) in ()HBj) . we 
obtain 



IGtt (27r)^ ™ \ r J J^i Jo cos^ lu y^{r,uj,m/ E) 



where we have defined 



ri^2 sin uj 1 
GTM^B) := / du; \_^ (54) 



A 
B 



j F{r) (55) 

l-(|)'i^(r). (56) 

§The elevation-angle of the radiation is defined as (7r/2 — tj). High elevation- angle means small zenith-angle 
Lo ^ and low elevation- angle means large zenith-angle uj ^ tt/2. 
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The definite integration in (j54j) can be performed as 



(57) 



wfiere 



a, 



■^\IB- Atan^u; 



(58) 



is a primitive function of the integrand. For convenience of the following calculations we show 
several evaluations of the primitive function ()58|) : 



Go 



'1 



-{ — ) 

F VrBH/ , 
1 / r \ ^ |r — 

0. 



F{r) 



r 



A 



my 
e) 



We rewrite the density-distribution Af{E,r; f) in (|H7|) as 

1 9 



-fT,,iE)AnE^ 



1 1 



4(27r)3^^«H^ ^ V r y F^{r) 

where we have defined a dimensionless positive function G as 



G 



m 



G(^,rj = Gi (^^,rj +Giia (^^,rj +Giib±,c (^,r) +Giib 



m 



(59) 

(60) 
(61) 

(62) 

(63) 



The functions Gi, Gna, Giib±,c and Gub are corresponding to the partial density-distributions 
A/i, A/iia, Mibit.c and A/iib in (jSZ|)- The relation between G's and A/''s is also given by (jH^ . The 
factors before a symbol of the first product in (jU^ correspond to the number- density of the 
particles radiated from the normal black body with temperature Tbh- The factor {r^u/f^y is 
a geometrical factor of spherical source of the radiation. The factor 1/F^{r) means a principal 
general relativistic effect and plays important role near the horizon. The factor G/2 means a 
residual correction including the finite mass effects. We can show that G/2 is regular function 
and have finite value (j(rBH)/2 = 1 on the horizon. We can also show that G is finite for all 
position r except for G{l,r — > oo). Therefore the density-distribution ()62j) is governed by the 
factor 1/F^(r) near the horizon. When m < E, the density-distribution (j62|) is governed by 
the factor {r^ii/ry for distant region. 

We concretely calculate the correction factors Gi, Gu^, Gnb±,c and Gub by using the rela- 
tions dMI), dSni) and dnH). From §^ we obtain 



rBH 



my ^ 

'ej 



- r 



(64) 
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for Tbh ^ r ^ oo. We find 



lim Gi(^,r 



tana;,2(f) 



dL 



and it is consistent with the relation 

hm 

between f)T7|) and From we obtain 



hm Voo-f^dE 



(m 



f) -1 







Wci <r <oo) 



(65) 



(66) 



(67) 



By using the property of the absolute- value function in and (jHTjl . we can summarize 



Ci+iia I — , T 



■) - «.g,0-«n.(|,.) 



r y 



'1 



my r - Tel 



1 



(68) 
(69) 



as the G-function for {] < m/ E < 1 and tbh < r < cxd. Both Gi and Gna as a function of r are 
non-different iable at r = rd, however, the function Gi+na becomes differentiable at r = rd. 
The function Gi+na is smooth and monotonously increasing function for tbh < r < oo. 
The function has the lower bound 2 for r = Tbh and the upper bound ()65p. The function 
G^i+iia is almost flat for r ^ tbh except for m/E = 1, then it can be approximated with 
Gi+iu{m/E,r) ~ Gi{m/E, oo) for r > tbh- 
From ()40|) we obtain 

(70) 



Ciib ( 



2 X 




1 - (f ) F(r) (rBH < r < (f 
(rc2 (f ) < r < oo 



This result has contained all contribution from the parameter region ^9/8 < m/E. This is a 
finite and continuous function of r. This is also differentiable except for r = rc2- 
From ()41|) we obtain 



Giibzt.c ( 



2 X 



/ r 



f ) ^(r) 



(rBH < r < r,! (f )) 



X < 



i_(=)V-t^J{ (1)^-1 



(t) 



< r < Tf, 



(f)) 



■ (71) 



l-(f) 



(res (f)<r<re2(f)) 



rc2 ^ < r < cx) 
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for 1 < m/E < y9/8. This is a finite and continuous function of r, however, is not differen- 
tiable at r = rd, res and rc2- 

When m/E = 1, the function Gi and G'iib±,c diverge for the hmit r ^ oo because the 
particle velocity at the infinite distance Voo becomes zero. However this divergence does not 
crucially affect the effective potential (fT7j) because the integration by E in ()18|) moderates the 
divergence and the effective Higgs mass (|18|) becomes finite for r — oo. 



6 Dynamical Formation of Spherical Domain Wall 

In this section we concretely derive the effective potential for the Higgs vev and consider the 
wall formation. The effective potential (fTTjl is governed by the effective Higgs mass /igg(r; A/") 
defined in (|TR|) . 



6.1 Effective Higgs Mass around Black Hole 

By substituting the particle-density-distribution M{E,r) resulted in (jU^ . the effective Higgs 
mass nl^ir-^N) in (fT^ becomes 

By defining a form function 

lift ' /'I ■> 

we can rewrite 



Hfi^^r] = — ds sMs)-G[^-,r], (73) 



/ 

The shapes of the form function ()73j) are shown in Figure IHl The form function ()74j] is defined 
for Tbh ^ r ^ oo. The range of the form function (ff^ becomes < Hf{m/TB}i,r) < 27/16 
for bosons and < Hf{m/T^ii,r) < 27/32 for fermions. 

Non General Relativistic Limit — If we omit the general relativistic (GR) effects as F{r) = 
1 and ^Gf = 1, we obtain Hf = 1 for any boson and Hf = | for any fermion. Then the 
effective Higgs mass becomes 

where we have defined the wall formation constant as 



/^eff(0 = V + — , (75) 



:= -^YgfYf (76) 



24 




Figure 6: Shapes of the form-function Hf{^^,r) defined in (fTHjl for various rrif/TBn and for 
(a) bosons and for (b) fermions. The horizontal axis d/r-BR is the distance from the horizon 
d = r — tbh, which is normahzed by the Schwarzschild radius tbh- For any ruf/TBR, we 
have the near-horizon hmits Hf{r —>■ tbh) —>■ 1 for bosons and Hf{r tbh) ^1/2 for 
fermions. For massless particles, we have Hf{r oo) 27/16 = 1.6875 for bosons and 
Hf{r —>■ oo) 27/32 for fermions. 
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and an effective degree of the freedom as 

{Of if: boson) , , 

29 f if ■ fermion) 

Tfie effective Higgs mass (f73jl reproduces the result in our previous work |H] where the GR 
effects including the red-shift-effect near the horizon and the gray-body-factor are omitted. 

High Temperature Limit — When we consider the high temperature limit Tbh rrif, the 
particles with high energy E ^ rrif dominate, then we can approximate 

\G{^,r) . lG(0.r) 

We can perform the ^'-integration in ()72|1 because the function G does not depend on E in 
the limit. We obtain Hf{0,r) = ^G{0,r) for bosons and if/(0,r) = jG{0,r) for fermions. 
The effective Higgs mass (f?^ becomes 

1 1 



/iefrW ^ -^i' + -G{0,r). (79) 

r [r) / 



We can approximate 



F{r) H y I (r ~ tbh) 

because |G'(0, r) is the slowly increasing function which takes value from 1 to y| = 1.6875. 
When we consider the distant region r ^ tbh, our high temperature limit ()80|) agrees with 
the non general relativistic limit ()75|) except for the factor The difference of the factor 
yI comes from the difference of the absorption cross sections, namely, the gray body factor 
makes this difference. When we omit the GR effects for the radiated particle, the absorption 
cross section is given by the horizon area Airr^^. On the other hand, the absorption cross 
section for the massless particle which obeys the Schwarzschild metric becomes ctbh = x^'^bh 
in pUj) . When we consider the near-horizon- region, the effective Higgs mass (jHUjl is governed 
by the factor 1/F{r) which diverges on the horizon. 

Low Temperature Limit — When we consider the black hole with the low Hawking tem- 
perature limit Tbh rrif, we can approximate 



Ig^'-G.. ^ j fe)V^-(^)'^W (rBH<r<r.(l^) 



(rc2 (^) < r < oo) 



2 2 

Then the effective Higgs mass (f7^ becomes 

1 1 



^ ^ f ^BH ' 



m'iF(r) , , 
- (82) 
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When we consider the distant region r ^ ^c2(5?^) where the particle with the energy E ^ Tbh 
cannot reach, the only high-energy- side of the tails of the distribution f^^^^ (E) contributes to 
()72|) . The difference between boson and fermion is negligible for ()82|) in the distant region 
because the high energy tails of the distribution fT^^{E) for bosons and for fermions are 
almost the same. Therefore we can adopt the Maxwell distribution to evaluate (|H^ in the 
distant region instead of the Fermi-Dirac or the Bose-Einstein distribution. By using the 
Maxwell distribution fxj^^{E) = exp(— £'/Tbh), we can perform the integration and obtain 



V + 



2 



Ki 



ruf 



BH 



BH 



3) 



87r2 F{rj ^ 

where Ki is the modified Bessel function of the second kind. On the other hand, when we 
consider the near horizon region tbh < r < rc2(^^), we can approximate 

1 1 



/^eff(^) 



A8F(r 



4) 



/ 



where we have defined the local effective degree of the freedom as 



9f{ 



g X exp 



-O.6OV2JF 



1 



-g X exp 



- 0.39\/2 



Tbh 



(r 



BH 



3/2' 



(boson) 
(fermion) 



(85) 



The only the radiated particle which can reach the position r contributes to the local effective 
degree of the freedom gf{r) on the position r. 

Summarized Form — Finally we can summarize the approximated effective Higgs mass 
for all parameter regions as 

A{rf 



1 

F(V) 



(86) 



where we have defined the local wall formation parameter A{r) as the same way of (j76|) by 
using the local effective degree of the freedom ^/(r) in (jH3j) instead of ^/ in (f77|) . As compared 
with the exact form ()74|) . the approximated form ()86|) contains an error of the factor one. By 
defining the local Hawking temperature as 

Tbh 



TBH(r) 



F{r) 



(87) 



we can approximate the local effective degree of the freedom ()85|) as follows 



9f[r} 



9f (TBR{r) > rrif] 
fTBH(r) < mj] 
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This property is very similar to the temperature-dependency of the degree of freedom in the 
relativistic thermodynamics. We present a energy-density-distribution p(r) near the horizon 
as p(r) ~ ■^g*f{TBR{r))T^ii{r) in Appendix B. This energy-density-distribution reproduces 
the relativistic thermodynamical relation with the temperature T = Tbh('")- The region near 
the horizon is not thermal equilibrium which is required for thermodynamical treatments, 
however, several behaviors near the horizon are very similar to that of the thermal equilibrium 
with the local Hawking temperature Tsnir). 



6.2 Wall Structure Solutions 



By substituting the resultant effective Higgs mass (IHUj). the effective potential (fTTj) becomes 



effl 



V + 



1 A\r) 
F{r) 



The equation for the Higgs vev in ()15|) with the effective potential 



becomes 



2 I + F{r) 



yZ 



(89) 



(90) 



The factor 1/F{r) in both the effective potential (jHn|l and the field equation ()9()|1 always 
diverges on the horizon. The form of the effective potential ()89|) near the horizon is different 
from that in the distant region (see Figure[7j). The minimum of the effective potential Ves{4>, r) 
in fl89|l depends on the position r. The effective potential Ves{4','r') is minimized by = 
near the horizon because the sign of /iefj(r) in ()8(ij) inverts when the position r approach the 
horizon. Therefore we expect the formation of the spherical wall-structure of the Higgs vev 
around the horizon even if Tbh < f^- It is difficult to define the exact radius of the wall- 
structure because the Higgs vev is continuously varying in the wall-structure. However we 
can define the characteristic wall radius row by 



fJ'esir Dw) — 0, 



(91) 



as a typical scale of the wall-radius. On the radius the sign of fils{r) inverts. 

To simplify the analysis of the wall structure we assume that the wall formation parameter 
is a constant. The characteristic wall radius becomes 



(92) 



as one of the solutions of the equation (jH^. The effective potential ()89j) is minimized at each 
point by the value of the Higgs field: 



V 





i^(^Dw) r^w 
F{r) r2 



1/2 



(r > tdw) 
(r < row) 



(93) 
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Figure 7: Function-forms of the effective Higgs potential Ves{(f),r) which is depending on 
the distance from the black hole. The dot in the potential-curve means the minimum of 
the potential at each r. The thick curve is the distribution of the Higgs field (pmmi'f') which 
minimizes the effective potential V^fj(0,r) at each r. 

The form of (/'mini'") is shown in FigurelHlfor various values of c/dw/'^bh- The normalized form 
of 4>min{r) is determined by dnw/i^BU- 

The physical phenomena outside of the horizon are considered and the inside of the horizon 
is not relevant to our subject. To clarify this we introduce a new positive coordinate d instead 
of r as a distance from the horizon: 



tbh + d. 



(94) 



We also define a characteristic- wall-distance ciow instead of the characteristic- wall-radius tbh, 
which satisfies 



(95) 



We introduce a profile function of the Higgs vev as 

1 



fis) 



VV2 



(rBH + d^ws) 



(96) 



which is dimensionless function and has value / = 1 in the ordinary vacuum. The parameter 
s is a distance from horizon normalized by (iDW- The equation of the Higgs vev ()90|) becomes 



d 



(so + sy ds 



A' 



2 (so + 1) 



(97) 
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Figure 8: Distributions of the Higgs field 0min('") which minimizes the effective potential 
Ves{4>,r) at each point. The horizontal axis is the distance from the horizon d normalized 
by the characteristic- wall-distance (icw- (We note the relations r = tbh + d and row = 
'"bh + c^DW-) The dotted line = is the ordinary Higgs vev and the vertical axis is 

normalized by f / a/2. The form of (pminir) is depending on the ratio (iDw/''^BH- The uppermost 
curve describes (pmini^) in the distant horizon limit d^w oo with keeping tbh constant. The 
lowest curve describes (pmini''^) in the near horizon limit c/dw ~* with keeping tbh constant. 
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where we have defined a ratio Sq := tbh/c^dw- This differential equation for the profile 
function is depending on both the ratio Sq and the wall- formation constant A^. 

To determine the Higgs vev structure around the black hole by the differential equation 
()97|1 . we should consider the boundary conditions. In the ordinary vacuum, the Higgs vev 
is = v/V2 which is given by the minimum of the bare Higgs potential (@]). The resultant 
effective potential Ves{(p,r) in (jHIJ|) agrees with the bare potential (jH) for r ^ oo. Then the 
boundary condition 

(a) 0(r) vj \/2 for r — > oo 

is required. This condition is equivalent to a condition f{s) 1 for s — >■ oo. If there is 
no Hawking radiation, we expect that nothing special about the Higgs vev happens on the 
horizon. When we turn on the Hawking radiation, we expect a finite deformation of the Higgs 
vev around and on the horizon. Therefore the boundary condition 

(b) 0(r) is finite on the horizon r = Tbh 

is also required. This is equivalent to the finiteness of the profile function f{s) on s = 0. 

Provided that the length-ratio Sq = tbh/ duw and the wall formation constant A'^ are given, 
the solution of the differential equation (jHIj), which satisfies the boundary conditions (a) and 
(b), is uniquely determined. In the Appendix C we consider the solutions near the horizon 
{d -C (iow) analytically. The near-horizon-solution which satisfy the boundary condition (b) 
is given by 



where C+ is a positive constant and 2F1 is the hypergeometric function. This solution is 
defined for s > 0. The near-horizon-solution ()98j] is a monochromatic increasing function and 
has property /(s) ^ for s — 0. We can numerically solve the equation (jHIj) so that the 
numerical solution is a numerical extrapolation of the near-horizon-solution (jUHj) . The value 
of the constant C+ should be so chosen that the numerical solution satisfies the boundary 
condition (a). Namely, the constant C+ is determined by the matching for the near-horizon- 
solution and the numerical solutions. The numerical solutions for various value of both the 
parameters r^n/dDW and are shown in Figure IHl 

The limit rfcw ^ '"bh in Figure M (a) is corresponding to the limit Tbh ^ because of 
both the relation d-ow — ^bw — A/ fi in this limit and the relation Tbh = In this limit 

the Hawking temperature Tbh is much greater than the critical temperature (~ fi) of the 
phase transition of the Higgs field. The numerical result shown in Figure IHl (a) indicates the 
formation of the wall-structure of the Higgs vev. The general relativistic effects are negligible 




(98) 
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(b) Comparable case of djy^ = ''^bh 
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(a) Wall with distant-horizon-limit 
1 




DW 

fc) Wall with near-horizon-limit 



Figure 9: Resultant structures of the Higgs vev around the black hole, (a) A wall with distant- 
horizon-limit ((iow ^ '"bh), (b) a comparable case c/dw = '"bh and (c) a wall with near- 
horizon-limit ((iow ^ '"bh) are shown. For each figure the solutions for = 0.001,0.01,0.1 
and 1 are shown as thin curves. The thick curves ^minl?") are the Higgs fields which minimize 
the effective potential at each point. 
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because the formed wall-structure is much distant from the horizon. This result reproduces 
the wall-structure concluded in our previous work jH] where the general relativistic effects are 
omitted. 

When we consider a comparable case d^w = '^bh in Figure|ni(b), the Hawking temperature 
Tbh and the critical temperature (~ ^) of the gauge- Higgs- Yukawa theory are also comparable. 
We also find the wall-structure around the black hole in the same way. The Higgs vev vanishes 
on the horizon: (pir^n) = 0, then the symmetry which has been spontaneously broken down 
by the bare Higgs potential is restored on the horizon. Therefore the wall separates the 
symmetric phase region on the horizon from the background of the broken phase vacuum 
(f) = v/ The restoration of the symmetry on the horizon is essentially caused by the 
general relativistic effect of the Hawking-radiated particles. 

The wall-structure for the near-horizon-limit c/dw ^ '"bh shown in Figure El (c) is more 
interesting. This limit is corresponding to the limit Tbh /W, namely, the Hawking tem- 
perature is much lower than the critical temperature (~ yu) of the phase transition of the 
gauge- Higgs- Yukawa theory. In the previous work [Hj where the general relativistic effects 
are omitted, we concluded that such a cold black hole cannot form the wall-structure by the 
Hawking radiation. However the general relativistic effects of the Hawking radiation form the 
wall-structure shown in Figure IHl (c) and restore the spontaneously-broken-symmetry on the 
horizon. The characteristic-wall-distance in this limit becomes 

^2 ^^ 

(iow ^ ^. (99) 

Finally, we conclude that the Hawking radiation of a black hole forms a spherical domain 
wall around the black hole even if the Hawking temperature is much smaller than the critical 
temperature of the phase transition. 



7 CONCLUSION AND DISCUSSIONS 

In this paper we have proposed a general relativistic (GR) formulation of the ballistic model 
to consider the Hawking radiation from the Schwarzschild black hole in the vacuum of the 
gauge- Higgs- Yukawa theory. We have found that the wall-structure of the Higgs scalar vev 
is formed even if the Hawking temperature is smaller than the energy scale of the gauge- 
Higgs- Yukawa theory. In the previous work where the GR effects were not considered, the 
critical Hawking-temperature for the wall- formation was found out |8j, however, we find out 
the absence of the critical Hawking-temperature. When the Hawking temperature Tbh is 
much smaller than the energy scale of the field theory ~ /x, we find the wall-structure closely 
near the horizon. The approximated width of the wall-structure near the horizon becomes 
the characteristic- wall-distance in ()99p. On the horizon we also find the restoration of the 
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symmetry which has been spontaneously broken down by the bare Higgs potential (0)). These 
results are essentially produced by the GR effects on the Hawking-radiated particles and on 
the Higgs field. 

Here we discuss validity of the GR formulation of the ballistic model. The ballistic model 
is valid when (a) the mean free path of the radiated particles is longer than the length- 
scale of the wall-structure dow and (b) the mean wavelength of the radiation is shorter 
than the length-scale (iow- The condition (a) is satisfied when the Hawking temperature 
Tbh is smaller than the critical Hawking temperature of the thermalization Tgy which is 
much greater than the critical temperature of the gauge- Higgs- Yukawa theory (~ /i). When 
Tbh is greater than Tgy, the wall-structure by the thermal phase transition rather than our 
mechanism is formed. This can be confirmed by considering the interaction-rates around the 
black hole jH E] . The mean wavelength of the radiated particles is approximately given by the 
Schwarzschild radius tbh for the region distant from the horizon, where the GR effects for the 
radiation are negligible. When we consider a situation Tbh ^ the wall-structure distant 

from the horizon ((iow ^ ^bh) is formed and the condition (b) is satisfied. On the other hand, 

when we consider a situation Tbh ^ A*, the wall-structure near the horizon (rfow ^ ''^bh) is 
formed and the GR effects become important. The mean wavelength near the horizon becomes 
short by the red-shift effect. Parikh and Wilczek showed that the Hawking radiation can be 
regard as the quantum-tunneling of the particle through the horizon [7^. They employed a 
particle description, which is very similar to our ballistic model, and the correct spectrum 
of the Hawking radiation is reproduced. In their argument the particle description near the 
horizon becomes valid because of the red-shift effect. Then we also expect that our ballistic 
description is valid not only for the distant-horizon but also for the near-horizon. 

In our GR formulation of the ballistic model, the Higgs scalar vev around the black 
hole is determined by the effective potential Ves{4>,r;J\f) in (fT^ which is depending on a 
density-distribution N'f{E,r) of both the particle-position r and the particle-energy E. The 
symmetry-restoration on the horizon and the wall formation near the horizon are essentially 
caused by the divergence of the density-distribution Aff{E,r) on the horizon. The density- 
distribution N'f{E,r) is depending on the particle-fiux ff{E,uj) on the horizon, which is the 
distribution of the initial radiation-zenith-angle u and of the particle-energy E (see Figure 
121). The divergence on the horizon of the distribution Mf{E, r) is caused by the cos~^ uj factor 
of the resultant flux ff{E,uj) in ()51|) . The flux ff{E,uj) on the horizon is quite different 
from that for the normal black-body-radiation ()52j) which includes the factor cosu. Therefore 
domination of the low elevation-angle (i.e., high initial radiation-zenith-angle uj ~ 7r/2) on 
the horizon plays important role in the wall-formation near the horizon. 

Here we summarize the derivation of the flux ()51|) with low elevation-angle domination. 
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The flux (|5ip is derived by an analytic continuation of the flux defined for high elevation- 
angle- region < uj < Uc into the flux for all elevation-angle < u < tt /2. The flux for the 
high elevation-angle < u; < c^c is calculated by the two assumptions; (i) the motion of the 
particles radiated from the horizon is determined by the geodesic on the Schwarzschild space- 
time and (ii) the observer at the infinite distance finds a disk-image with uniform intensity of 
the black-body-radiation for each particle-energy as the image of the radiated particles. The 
assumption (i) is the simplest application of the general relativity to the Hawking-radiated 
particles. The assumption (ii) may be valid because the observation of the image with non- 
uniform strength means that the Hawking radiation is not thermal radiation and carries out 
some information from the black hole. 

The domination of the radiation with low elevation-angle is consistent with the picture of 



the heat-bath around the horizon with the local Hawking temperature TBu^r) ~ Tbk/ y F{r) 
(see Appendix A). The energy-density-distribution of the Hawking radiation p(r), which is 
derived from the resultant particle-density-distribution Aff{E,r) in (jHH), is also consistent 
with TBuir) (see Appendix B). We have pointed out that Tbh(^) naturally arises in the 
resultant effective Higgs mass in (jHSj). Hotta discussed that the temperature closely near the 
horizon becomes very high due to the picture of the heat-bath with the local temperature 
2^bh('") and the thermal phase transition in the string theory arises |14j . 

The total energy of the Hawking-radiated particle between the sphere near the horizon 
with a radius tbh + d and the distant sphere with a radius R is given by 



This total energy should be finite, however, this has an IR divergence for R ^ oo and an 
UV divergence for d ^ 0. Because the black hole has a finite lifetime tbh, there arises a 
natural IR cut-off for the radius -Rcutofr ~ tbh Then the IR divergence is easily solved. On 
the other hand, the UV divergence -Etotai ~ l/rf is an open problem. A cut-off of the Hawking 
radiation near the horizon may be required. This cut-off is equivalent to the cut-off for low 
elevation- angle domination, namely, the analytic continuation is restricted toO<Ci;<7r/2 — e 
with finite e. This problem may solved by quantum gravity, then the cut-off may be given 
by the Planck length scale dcutos ~ ^pi- If we take the cut-off in the Planck length scale, our 
argument in this paper is valid for d > l^i. Therefore we expect that our results are valid as 
long as we consider low energy (^ mpi) phase transition. Hotta discussed that the cut-off is 
given by Tbh ('"cutoff) ~ "^pi and all of energy and all of the information (entropy) are carried 
by the Hawking radiation, i.e., -Etotai = "^bh [El- In this case, the cut-off rcutofr is much 
smaller than the Planck length Zpi. 

Finally we discuss applications of the wall-structure. We expect a net charge-transportation 
into the black hole by the Hawking radiation and a phenomenon of the spontaneous charge-up 





(100) 
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of the black hole by the effect of the wall-structure |HI • We also expect the baryon number pro- 
duction [TJ by the transportation of the hyper charge ^21 El- This production-mechanism 
may be realized by the thin-wall black-hole baryogenesis or the direct black-hole baryogenesis 
proposed in |n|. 

The wall-formation near the horizon is caused by the Hawking-radiated particles with low 
elevation- angle. All of such particles return into the horizon and are not directly observed by 
the observer at the infinity distance. One may consider that the radiation with low elevation- 
angle and the wall-formation near the horizon belong to a metaphysical subject. However 
the formation of the wall-structure near the horizon can be confirmed by the observation of 
the net charge flux from the black hole. Therefore the wall-formation near the horizon is not 
metaphysical. 
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APPENDIX A: Radiation Angle Distribution and Met- 
ric 

In this appendix we show that the radiation-angle-dependency of the the Hawking radiation 
on the horizon (the factor 1/ cos'^ uj in (jKTjl ) is naturally derived from the property of the space- 
time. Around the black hole described by the Schwarzschild metric a short line-element 
dr near the horizon in the Schwarzschild coordinate system becomes a long line-element 

dsr = -^dr. (101) 



in the proper coordinate system. Namely the Schwarzschild coordinate system compress the 
r-direction near the horizon (see Figure ITUI) . This compression changes the zenith-angle near 
the horizon. The zenith-angle of a line-element in the the Schwarzschild coordinate system 
(the thick segment in Figure ITUI fB)) is defined as 

rdO 

tana; = — (102) 
dr 

which is used for the radiation-zenith-angle in the main part of this paper. On the other 
hand, the zenith-angle of the same line-element in the proper coordinate system (the thick 
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segment in Figure ITUl (A)) becomes 



rd6 

tan(^ = — . (103 

dSr 



These zenith- angles u and <f are different description of the same hne-element, then a relation 
between these angles is derived from the definitions (jl(J2|) . (jl(J3p and the relation (jl(Jl|) as 



tamp = vFtana;. (104) 
For convenience we define a function 

A{F,u) -.= 1+ (^j-l^cos^u, (105) 



then we obtain several relations 

d(j) 1 1 



(106) 



9 , 1 cos^o; , _, 

9 , sin^ u , , 

We put a angle-distribution of particle-flux on the (/^-coordinate system as 

dT = f{E,(p) X dE X dipsiiKfi dtp. (109) 
The angle-distribution ()109j) can be rewritten for the tu-coordinate system as 

dJ^ = f{E,(p{uj))^ ^^^'^^^^ X dE X du sin 00 dip 
did sin uj 

= f{E,uj) X dE X dusinu} dip, (HO) 

then we obtain the relation 

f{E,u) = fiEM-)) py,^lp^^^,/, . (Ill) 

which gives us the coordinate transformation law (u; ^ if) for the angle-distribution. 

When we consider the near-horizon-limit r tbh, "we have F{r) and the coordinate 
transformation for the zenith-angles (jl04p becomes 

, . Jo (0<a; < 7r/2) 

^ { ■ (112) 

[ 7r/2 {u = 7r/2) 

We also obtain the following forms in the limit: 

A{F,u) —cos^u, (113) 
F 

f{E,cu) ^ f(E,^{u;))^^. (114) 

COS'^ UJ 
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(A) Proper Length Element and Angle (B) Coordinate Length Element and Angle 



Figure 10: (A) A line-element in the proper-length coordinate system (thick line segment) and 
(B) the corresponding line-element in the Schwarzschild coordinate system near the horizon 
(thick line segment). The relation between the coordinate line-elements is dsr = F~^^'^dr. The 
definitions of the zenith-angles of the line-element are tamp = rdO/dsr and tancj = rdO/dr. 
Then we have the relation iamp = a/F tana;. 



If we assume that the zenith-angle-distribution of the fiux /(-E, 9?) in the proper coordinate 
system is smooth around = 0, we have 



f{E,u;) = f{E,0) 



COS^ 00 



(115) 



in the near- horizon-limit. Here we find another derivation of the factor l/cos'^tu in the 
angle-distribution of the particle-fiux (jllSj) . In Section |3] we have derived the distribution of 
the particle-fiux of the Hawking radiation in (jHTjl including the factor 1 / cos'^ u by analytic 
continuation and uniformity of the Hawking radiation. On the other hand, the factor 1 / cos^ lo 
is derived from the smoothness of the distribution and geometry of the black hole in this 
appendix. 

Inversely if we assume the resultant fiux (|5ip. we obtain the fiux on the proper coordinate 
system as 

1 1 



(116) 



where we have defined the proper energy of the particle E{r) := E/J2F{r) and the effective 



local Hawking temperature TBu{r) := Tbh/ j2F{r). The proper fiux for the zenith-direction 
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~ resulted in is consistent with the heat-bath picture with the local Hawking 

temperature T^nir). 



APPENDIX B: Energy Density Distribution around the 
Black Hole 

In this appendix we calculate the energy density of the Hawking-radiated particles around 
the black hole to clarify the physical meaning of our resultant particle number-distribution 
Nf{E,r) in (jU^ . The energy density of the Hawking-radiated particle at the position r with 
the particle-species / is given by 



Pf{r) = / dEEAff{E,r). (117) 



By defining the form function 



^BH 



we obtain 

P/W - TS«/^B„x^^^j A-^g^.^j. (119) 

We display the shapes of the form function Kf{r) in Figure ITTl For any mf/T^n, we have 
the same near-horizon limits Kf{r — > tbh) — * 1 for bosons and Kf{r —* tbh) — > 7/8 for 
fermions. For the massless particle we have Kf{r —>■ oo) ^ 11 — -'-•69 i^i bosons and 
Kf{r — s> oo) ^ II X I ~ 1.48 for fermions. 

The energy-density flll9p has an interesting approximated form: 

2 2 

P/W ^ l^9*fiTBn{r)) fUr) x (^) , (120) 
where we have defined the local Hawking temperature 

TBH(r) := ^^f= (121) 

and also defined the degree of freedom of the particle / 

9fiT) if ■■ boson 
l9fiT) if : fermion) 

for the temperature T with the fermion- correct ion. The relation ()120j) is just the thermody- 
namical relations between the temperature Tbh('") and the energy density p/(r) except for 
the spherical geometrical factor ir-Bn/r)'^- Then our evaluation of the differential particle- 
number-density fl62|) is consistent with the picture of a heat-bath around a black hole with 
the local temperature T(r). 



9*fiT) ■■= { r/ \ (122) 
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Figure 11: The shapes of the form-function Kf{^^,r) defined in ()118|) for various irif/T-Qn 
and for (a) bosons and (b) fermions. The horizontal axis d is the distance from the horizon 
of the black hole with the Schwarzschild radius tbh- 

APPENDIX C: Solution near Horizon 



We consider the properties of the solution closely near the horizon: r = r-QYi + d with d <^ tbh- 
In the region near the horizon the positive contribution to the Higgs mass term dominates 
in the effective potential (jH^ because the Schwarzschild factor becomes very large near 
the horizon^. The field equation (jHH) for the profile function /(s) = :;;/^0(^bh + c^dw*) is 
approximated to the linear differential equation: 



d 



(so + s) ds 



[so + s)s -^f 



~2~s 



/, 



(123) 



where we had defined the ratio sq 
differential equation is given by 



'"bh/c^dw- The general solution of this second order 



f{s) 
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(124) 



^The reason can be regard that the local Hawking temperature T^yl/ \/F{r) becomes very large. 
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where 2F1 is the hypergeometric function. The coefficients C+ and C_ are determined by the 
boundary conditions. The ffist term is a monochromatic increasing function and has zero 
value on the horizon {s = 0). The second term has singularity on the horizon. To obtain 
the Higgs vev by the differential equation ()12H|1 we have required the boundary conditions; 

(a) f{s) 1 for s 00 and (b) f{s) is finite for s —* 0. Due to the boundary condition 

(b) the second coefficient C_ should be zero. The ffist coefficient is determined by the 
matching of the solution ()124|) with C_ = and the numerical solution which satisfies the 
boundary condition (a). As a result of matching to the numerical solutions shown in Figure 
ini we obtain 0(1) values of C+. For example we obtain C+ = 0.90 for = 0.001 with the 
wall extremely near the horizon, Cj^ = 1.00 for = 0.001 with c/dw = ''"bh and C+ = 1.14 
for = 0.001 with larger wall than the scale of the horizon. 

The analysis of the Higgs vev near the horizon results that the symmetry broken-down 
spontaneously by the bare Higgs potential is restored on the horizon /(O) = and that the 
spherical wall-structure of the Higgs vev separates the symmetric phase region on the horizon 
from the broken phase vacuum in the region distant from the horizon. Therefore the resultant 
wall-structure is just the "domain wall" , according to the ordinary terminology. 

To clarify the behavior of the solution p24j) . we will consider the extreme cases. When we 
consider the wall with the near horizon limit, i.e., sq = tbh/c^dw 00, the solution becomes 
a simple form: 



When we consider the wall which is extremely distant from the horizon, i.e., Sq = tbh/c^dw 
0, the solution becomes 



In both cases, C_ = is required by the boundary condition (b). When the wall-formation- 
constant is much small, i.e., -C 1, the solution near the horizon which satisfies the boundary 
condition (b) can be summarized as 



When we consider the electroweak wall-structure, we have the wall-formation-constant ~ 
0.001 [S] then we can employ this form. 
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